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ABSTRACT 


Infrared vibrational-rotational transitions in a homonuclear 
diatomic molecule, can only occur at high pressures and/or in the presence 
of an external electric field. 

Experimental attempts at making a high pressure H, laser by using 
intense 8 particle radioactive isotopes as a source of preionizing elect- 
rons are described. 

Resonant Electron-to-Vibrational transfer is also discussed asa 
possibility for creating a population inversion in H,- An experiment 
designed to achieve E-V transfer, from metastable NO, electronic states 
to Hy vibrations, is described. Experimental values of the self quenching 
cross-section together with a velnetor quenching of NO, by H. are presented. 

Finally, non-resonant Electronic-to-Vibrational energy transfer is 
discussed. It is shown that resonance is not an important criterion 
for efficient E-V transfer. The percentage of vibrational energy produced 
in the quenching of Na! by H, is calculated, and compared with experimental 
values available. Non-resonant E-V transfer theory is discussed in light 
of the Born-Oppenheimer Approximation and it is shown how the latter offers 


a broad basis for selection of suitable pumping atoms or molecules for any 


potential laser molecule, including H,- 
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CHAPTER I 


INTRODUCTION 


High power tunable infra-red lasers would find many uses in 
research and industry. For example, they would be particularly useful 
in molecular spectroscopy and biology [1,2], in the field of communications 
[3] and in isotope enrichment applications [4,8,9]. 

Output pulses of very short duration (of the order of picoseconds) 
could be obtained by mode locking such lasers. This follows from the 
fact that in mode locked operation the output pulse width is inversely 
proportional to the bandwidth [16], which in a continuously tunable 
range,could be very large. 

Many tunable semiconductor lasers have successfully been put into 
operation. These presently cover the region of the IR spectrum 
between about 0.5um to 50um [3]. However, the power output levels from 
such lasers are very low (not exceeding a few dozen watts at the most). 
This limits their uses considerably [1]. 

On the other hand, high output power IR gas lasers have been built 
in recent years. Output powers of up to a few Kilowatts continuous and 
a few Megawatts pulsed are relatively common from the co, laser. This 
makes tunable IR gas lasers a particularly attractive proposition. 

Basov et al [1,5] have pointed out the possibility of making tunable 
CO and Hy lasers by utilising pressure broadening of the spectral lines 
so that adjacent lines overlap, producing a continuous amplification band. 

Continuously tunable high pressure CO. lasers have been reported 
by T.Y. Chang and 0.R. Wood [10], by A.J. Alcock et al [11] and by 


N.G. Basov et aZ [12]. Chang and Wood used the 4.2ym radiation from a 
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Hydrogen Bromide pulsed laser to optically pumped pure CO, at 33 


Z 
Atmospheres. Alcock et al used UV photopreionization to obtain continuous 
tuning over a range of about Deena at 15 Atmospheres. Basov et al 
obtained operation around 50 Atmospheres by using the e-beam technique 
resulting in a tunable range of about LOGens 

For a CO laser a pressure of between 10 to 15 Atmospheres would 
be high enough to produce a continuous band [1]. 

Legis shown below that a pressure of about 40 Atmospheres would 
suffice for the H, laser. A tunable Hy laser has not been reported yet. 

A Hydrogen laser operating on vibrational rotational transitions 
would emit radiation in the near IR (at around 2um) [18]. 

(eacnulan et al [6] have achieved lasing action in H, and Dy around 
the lum region. Population inversion between the vibrational rotational 
levels of the first two excited electronic states was obtained by making 
use of the Franck-Condon Principle and the fact that the minima of these 
two states occur at different internuclear separations. The method of 
excitation was by an electric discharge at 35KV and 20Hz at relatively 
low pressures [7]). 

Some IR lasers together with their output wavelengths are listed in 
Appendix I. 

IR Hy laser transitions between vibrational-rotational levels would 
not occur under ordinary circumstances because, being a homonuclear 
diatomic molecule, Hy does not have a permanent dipole moment, and the 
dipole transition probability is zero, such transitions being "forbidden" 
[13]. An electric dipole moment could, however, be induced in the molecule 
by an external electric field (d.c. or a.c.) [5] or by intermolecular 


collisions [14] giving a nonzero, although small, probability for such 


transitions. Basov et al [5] have suggested that the strong electric 
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field of CO, laser radiation combined with the isotope effect in the 
unsymmetrical HD molecule (which has a dipole moment in the non zero 
vibrational quantum states [13]) could be used to make a high pressure 
Hydrogen laser. Christensen and Greenfield [14] have analysed a 
hypothetical H, laser that utilises both the collision induced dipole 


(cid) and the electric field induced dipole (efi) effects. 


1.1 Scope of this Work 

A condition for continuous tunability in a gas laser is derived 
below. It is shown that a Hydrogen laser operating at a pressure around 
40 Atmospheres would be continuously tunable. 

Chapter II deals with experimental attempts at obtaining a high 
pressure discharge in Hydrogen and other gases by using an intense 
radioactive 8 particle emitter for preionization. 

Chapter III describes an experiment designed to optically pump 
the Hydrogen molecule by resonant E-V transfer from NO, excited metastable 
states. Experimental cross-section values for the quenching of these 
NO, states by collision with other NO., molecules and with H molecules 
are presented. 

The general theory of E-V transfer is described and applied to some systems 

containing Hydrogen molecules in Chapter IV. ‘The results obtained in 
Chapter III are discussed in light of this theory and it is shown that 
resonance between electronic and vibrational levels is not as important 
a condition for efficient E-V transfer as is generally thought. Cross-- 


sections for the E-V transfer process are calculated for a few systems. 


These calculations are used as a basis for suggesting some atoms that 
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could transfer electronic energy to Hydrogen vibrational levels. 
Finally a summary is presented in Chapter V together with some 


suggestions on how the present work might be extended. 


1.2 Condition for Continuous Tunability 


When a radiating atom or molecule suffers a collision with 
another particle, the perturbation causes an interruption in the emitted 
“wave train". A Fourier transform of this train reveals additional 
frequencies that would not be present if no collision had occurred [17]. 
The emitted line is therefore ("collision" or "pressure") broadened. 
If the consequent increase in line widths is large enough, adjacent 
lines overlap and the emitted spectrum would then become a continuous 
baad. 

It can be shown that the profile of a pressure broadened line 


is Lorentzian [17] and is given by 


ue 


5 (1.1) 
(Av) -b ¢ 


g(v) = c¢ 


Pak ) 
fe) 


where 
c is a constant incorporating the oscillator strength of the 
particular transition under consideration, 
v is the frequency, 
vis the centre frequency of the profile, 


T is the mean free (or "flight'') time between collisions 


fo) 
and 
Av = v-v 
fo) 
(g(v) is symmetrical about and has a maximum at Vy)° 
From (1.1) it follows that the half intensity linewidth Av? is 
given by 


ee ane (1.2) 
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Also, from the definition of Ty? it follows that the average 
number of collisions per unit time, denoted by S, is given by 


ie (1.3) 


i 
ak 
° 

If the cross-section for broadening collisions is oO» v the 


mean relative speed of the particles and N their concentration, then 
vN Cha) 


From (1.2) and (1.4), it is seen that the half line width is 
proportional to the gas pressure at a given temperature (This assumes 
that oF is not a strong function of pressure). Further, in a gaseous 
mixture the half width of a line is equal to the sum of half widths due 
to each component of the mixture, considered separately. 


Equations (1.2), (1.3) and (1.4) can be combined with the ideal 


gas relation [19]: 


P = NRT (755) 
to give 
Av 
hw iL 1 
= eae —_ 1.6 
P Fe PE oe) 


For a given temperature, this ratio is approximately constant ina 


given pressure range, if o is regarded as the mean collision cross-section 


in this range. 
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For a gas like CO the comstant ratio in (1.6) is of the order of 
pote core [1]. 

Now v is inversely proportional to the paare root of the mass 
of the colliding molecules; this follows from the kinetic theory of 


gases [17]: 


gxr . 1/2 


v= ( ™ Gp) 


We deduce, therefore that for Hy the broadening per unit frequency, 
M 1/2 
(Av, Ps would be approximately 10! (CO ) where M__ and M, are the 
2 


M, co 


molecular weights of CO and H neepeceieaiy. This yields 


Z 


AV 


P 


a307 x 10’ Hz/torr (1:28) 


for H,- 

We will use equation (1.8) to estimate the pressure at which 
continuous tunability may be achieved in a Hydrogen laser. It is first shown below 
that to achieve overlap between the adjacent vibrational-rotational bands, 
it is only necessary to broaden the lines by an amount approximately 
equal to the rotational constant, B, of the molecule. 

A diatomic may be regarded as an anharmonic oscillator, the anharmonicity 
being due to the deviation of the potential curve of the molecule from 
the simple harmonic potential. This results in an uneven spacing of the 
vibrational levels, the higher ones being closer together than the lower 
ones [20]. In the Morse-Oscillator model (where the potential curve 


under the influence of which the diatomic is vibrating is assumed to be 


the Morse potential) the magnitude of this anharmonicity AE, is given by [1] 
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Neen pee (1.9) 
where 


w is the (radian)frequency of the 1+0 transition, 


0 
D is the dissociation energy of the molecule. 
The lines due to adjacent vibrational transitions do not, ordinarily, 


coincide because of this anharmonic defect. The separation (in 


frequency) Av, between them is given by 


If each line in this wibcationel spectrum were broadened by this amount, 
the adjacent vibrational levels would overlap. The anharmonic defect 
is of such a magnitude that it would be necessary to have pressures of 
the order of hundreds of atmospheres for the desired overlap to occur. 

Fortunately however, since each vibrational level has a rotational 
fine structure, in practice it is not necessary (for the case of H, and 
CO at least) to broaden the lines by AE. 

Because of the statistical profile of the kinetic energy distribution 
amongst the molecules of a gas, it turns out that each vibrational level 
is accompanied by a spectrum of rotational levels excited in the frequency 
interval AQ = “e [1] i.e. in the energy interval AE, = KT. At 300°K, 
this works out to be about 3 x 10°" CV Le AED " AE then we may expect 
overlap between adjacent vibrational levels even if each rotational 


line was broadened enough to overlap only the next rotational line in the 


same vibrational level [1]. For Hy and CO, the value for the anharmonic 
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defects are found from (1.9), to be about 5 x 10° ev and 2 x 10° ev. 


(The values for w and D being obtained from Herzberg [18]). 


LO 
The value for the anharmonic defect in Hydrogen (5 x Tom ev) 
is of the same order of magnitude as the value for spread in excited 
rotational energy levels (3 x none ev) (see also Appendix II), The value 
for the defect in CO (2 x 107 ev) is less than AE, (Sex ‘iy ev). 
Hence, for both these molecules it is only necessary to broaden each 
line by a value approximately equal to B, the rotational constant, (this 
being, approximately, the separation between rotational lines). 
The values of B for Hy and CO are 3.6 x tor? ev and 2.4 x tom ev 
corresponding to the frequency spreads of to" Hz and 10." Hz respectively. 
Finally, using the values 4 x 107 Hz/torr (from 1.8) and 107 Hz/torr 
(from Basov et al [1]) for these two molecules we find that neascites of 


about 40 Atmospheres and 13 Atmospheres, respectively would suffice to 


produce continuous vibrational rotational spectra in these two molecules. 


1.3 Condition for Population Inversion 


Homonuclear diatomics like Hy have long vibrational-translational 
relaxation times because dipole transitions are "forbidden". Rotational- 
translational relaxation times are, at 300°K and up to several atmospheres, 
many orders of magnitude smaller. Hence, when considering the inter- 
action of these three types of energies (translational, vibrational and 
rotational) it might be assumed that the (translational) temperature of 
the aera equal to the rotational temperature T. (as defined below). 


[T may be defined in the kinetic theory of gases by equation (1.7) above] 


[5]. 
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The number of molecules in the rotational level J, in the lowest 


vibrational state is given by [18] 


NhcB BJ(J +1) 
KT, (23 +1) exp{ Te, 


N@I) ‘= (1.10) 
where 
N is the total number of species: c the velocity of light and the 
other terms have been defined Peters’ 
This Boltzman distribution defines the rotational temperature, Te 
Further it is generally a good approximation to equate z witht [5]; 


The vibrational temperature, 7 is defined as follows [5]: 


ane) eee AE /RT,) | (1.11) 


N 
at 
where 
Ny and N) are the concentration of molecules in the first 
vibrationally excited and the vibrational ground states 
respectively, 
AE is the energy separation between these two levels. 
Now, vibrational-rotational transitions in molecules which have 
had a dipole moment induced by either the efi or cid effects must satisfy 
the following selection rules: [14] 
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where 


V and J are the vibrational and rotational quantum numbers, 
respectively. (The Raman transitions follow the same 
selection rules), 

Consider such a transition between two devele J, and J, such that 


2 1 


I = 4,42 (1.14) 


N(J,) Ny, (23,41) -B 
NOW Le NS Cue) expl gp tJ, 541) - J, (J) 4+1) 3) 
1 I 1 r 
N, (25,41) AE. 
Sa) ore er CRD | aera 
Ny (25,41) Kr 


where 
AE. is the energy difference between the two rotational levels. 
(We have assumed that the rotational constant B for both the vibrational 


N 
2 , 
states is the same). Substituting for yp from (1.11) gives 


N(J,) AE AE AE 

2D Vv 1 Vv r 

= exp( - =~ ) exp = [- =~ - =— (1.15) 
N(J,) KT K as ae 
or 

N(J,) AE 
pee ae aa ee (1.416) 
NG,) ~ exp{ KI se eACES) | 
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(25,41) 
where we have used (1.14) and the approximation ————~ = 1 which holds 


(25, +1). 
for the higher rotational levels. From (1.15) and (1.16) we’ can write 
the condition for population inversion [N(J,) > N(J,)] in the following 


equivalent forms [5, 14]: 


AE. 
or 
qT, AE. 
7 2 C18) 


ee the approximation TL = T has been used. 

Equation (1.18) implies that no Q branch (AJ = 0) transitions could 
ae because for such transitions AE. = 0. 

However, if in the transition process translational energy plays 
a part, as is the case for cid transitions, then equation (1.18) is 


modified as follows [14] 


i AE 
T 2 Gey DE) oe) 
( r t 
where AE. is that change in translational energy that accompanies the 
particular vibrational rotational transition under consideration. 
This means that a population inversion may exist even when AE = 0. 
Further, this inversion could be significant at a low temperature, T. 


This “translational partial inversion" allows Q branch transitions [14]. 
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It should be pointed out that conditions (1.18) and (1.19) have 
been derived assuming that the vibrational levels are mopulaved according 
to the Boltzman distribution (1.11) This is true if the molecules are 
considered as harmonic oscillators. In that case,Landau and Teller have 
shown that vibrational relaxation of the molecules may be described 
by one overall relaxation rate [46]. Hence, although the relaxation times 
of different vibrational levels of the harmonic oscillator are different, 
they fit together in such a way that one overall rate relaxation time 
could be defined for all the transitions. 

Treanor et al [64] have shown that if anharmonic effects are 
considered, the distributions could show substantial departure from the 
Boltzman scheme. Under these conditions, the Landau-Teller model is not 
quite accurate. 

Fortunately, it turns out that the anharmonic effect is advantageous 
to the process of establishing a population inversion especially in the 
important case of low translational temperature, since then the populations 
of the lower vibrational states can be erreiderahie smaller than those 


predicted by the Landau-Teller model [64]. 
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CHAPTER II 


ELECTRIC DISCHARGE EXPERIMENTS 


Pumping a laser gas by creating a discharge has been demonstrated 
to be an efficient technique. High powered efficient co, lasers have 
been successfully put into operation by this means. Seguin and Tulip [21] 
obtained increased discharge stability and large discharge energies by 
incorporating U-V spark radiation for initiating and sustaining the 
discharge. This method was later used by Alcock et al to obtain lasing 
action at pressures up to 15 Atmospheres in co, [11]., Tulip, Seguin and 
Faszer [22] used the technique to obtain a high repetition rate (up to 
4000 pulses per second) TEA laser. 

An extensive quantitative investigation of the radiation from a 


spark gap and its effect on the CO, laser mix was carried out by Seguin, 


2 
Tulip and McKen [23]. It was found,for example, that the absorption 
length of the U-V radiation at 200 torr (in a 1:1:1 and a 1:1:3 ratio 

with tripropylamine additive) was about 5 cm. With some other additives 
it could be as long as 35 cm. 

If the preionization technique is to work at high pressures (at 
several tens of Atmospheres) then it would be necessary to use much harder 
radiation to obtain the required longer propagation lengths. It can 
be shown for example, that soft x-rays in the region between 0.5A and 
UG would have propagation lengths between about 30 cm to about 2 cm at 
50 Atmospheres in air. Such wavelengths could be obtained from an x-ray 


tube operating between 20 and 50 KV. (X-rays produced below 70 KV 


dissipate most of their energy by ionizing the absorbing molecules [24]) 
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Another source of electrons for preionization would be a radio- 
active 8 particle emitter like Strontium 90 (er). The daughter product 
of this isotope, Yttrium 90 Ce) is also a 8 particle emitter. The 
average 8 particle energy from Strontium 90 is about 0.18 MeV. Such 
electrons would have an average absorption length of about 2 cm in 50 
Atmospheres of air [26]. The 8 particles from ie are on the average 
about four times as energetic as those from wart and would therefore have 
an absorption length of about 8 cm under the same conditions. The 
absorption lengths in any laser mix would be different from these figures. 
However, in the absence of specific data for such gases as H, and CO,» 
the absorption lengths for air, mentioned above, were used as the basis 
for selecting gol for the experiments described below. The other main 


reason for selecting this particular isotope was that it could be 


purchased commercially. 


Belk Were Experiment 


As pointed out previously, a Hydrogen infra-red laser would have to 
be operated at a bias EF) high enough for the collision induced dipole 
moment to be substantial. The pressure chamber (shown below) was designed 
to withstand pressures up to 50 Atmospheres. 

At these high pressures (tens of Atmospheres) a high voltage pulse 
of small time duration and a strong source of preionization would be 


necessary. Each of these factors is dealt with below. 


Pe elleeal Voltage Requirements 


The magnitude of the electric field strength required was calculated 
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using the Raether Breakdown Criterion [28] for Helium and air, data on 


Hydrogen being insufficient. See Figure 22). 


The criterion can be written as 
a d. = 17 + ind | QoL) 


where 
a is Townsend's first ionization coefficient (in Ey 
d is the (critical) avalanche length (in cm) at which the 
avalanche space charge field becomes comparable to the 
externally applied field. 
The numerical value given on the right hand side of equation (2.1) 
varies slightly with different gases [28], but the variation is small 


enough not to seriously affect the estimated value for the field. 


Equation (2.1) was used in conjunction with the Townsend equation 


= A exp (-Bp/E) (2.2) 


~|2 


where 

p and E are the pressure and field intensity respectively, 

A and B are constants depending on the gas and the value of E/p. 

The breakdown field intensity for He at 50 Atmospheres, using 
equations (2.1) and (2.2) is about 120 KV/cm. This figure was used as 
a rough target for the Marx generator that was made for the experiment 
(see Appendix 4). The latter was successfully operated to give output 


voltages higher than 150 KV. 
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The voltage pulse duration is important in any gas discharge. 
Rogoff [30] has shown that thermal instabilities develop in a 500 torr 
Hydrogen discharge in a time interval of about 100 ns leading to arcing. 
In the experiment an attempt was made to decrease the pulse width as 


much as possible by reducing the inductance loop areas. 


2.1.2 Preionization Requirements 
An approximate condition for the minimum preionization density, no» 


required to initiate a glow discharge is that the avalanche head's 


=1/3 


radial extension, r, be greater or equal to the dimension (n,) [sy 


ras (an? (2.3) 


where 
dq. is the extent of the avalanche along the field direction at 
breakdown and 4 is the mean free path of the electrons [32]. 
Hence the criterion for minimum preionization density is 


lyfe 


~ , \1/2 
(n,) < Q 4,) (2.4) 


Since } is inversely proportional to pressure this may be rewritten 


as 


-1/3 a 17.2 
Cpe it ) (2.5) 


(n,) 
where 
2X, is the mean free path at atmospheric pressure and (n,) the 
r Fr 


A 
minimum preionization density at P atmospheres. [>— has the 


unit of length since P is a number in (25) Ve 
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For a co, TEA laser Palmer [31] finds a, to be about 10° nee 
Hence at 50 Atmospheres ny =~4x 10° ene 

The ele sources used each had a dosage of 80 mRad (ceeie. (By 
definition, 1 Rad is the unit of dose that corresponds to an absorption 
of 100 ergs in one gram of absorbing medium [33]). Each source irradiated 
about 50 ene in the discharge chamber. Using the value of 15.4 ev as 
the ionization energy of Hy» it follows that the 8 sources produce 
electrons at the rates of 10° (Carey aec) and 10! (ane (sec) in 
Hydrogen at pressures of 50 Atmospheres and 1 Atmosphere respectively. 

These figures ignore secondary ionization processes like emission 
from the cathode. Losses due to attachment, on the other hand, are not 
considered either. 

The steady state electron concentrations would be close to these 
numerical values (i.e. 10° ae and 10! nat) since the 8 radiation from 
the isotopes is continuous. 


It is seen therefore that the gee source would probably provide 


enough preionization density for a discharge in Hydrogen at 50 Atmospheres. 


2.2 Experimental Procedure 


The experimental set up is shown in Figure 2.2 below. The Marx 


Generator is described in Appendix 4. 


The details of the Pressure Chamber are shown below in Figure 2.3. 


Sen pas amy ‘eee 


Pu 


hn © 


eo heel) aig GE Ba ea iy 
reiing acta fap a eaberrbley (ANY aoa te ‘bbe, 
berehfasvk sinew dived? TEC) wre acy ate rete ty. 
ne ge #08 Yo welch age goied .sodieiule eral 7 ; 
ST savboty @iseime A Aes cele sePet a @ ogee 
“ae “Goma” fe Tee aun coo aay Boy ae 
ai fea otacenaia S be no eee) a one 
rebunien mth sha wasarg peak Lack: qtahaasd: # A 
we oth ~Renl td70 aay uv  doaastion 90! ota , 


ewig 64 Geeto © rar, one soe ram bps 
want neweniaay € v4? ani ‘ ‘o? Soke’ ial 
ry 


ebay italy A ere Worn a yaa ips a i- 
re sugaianth O61 32 equ ihe! at < jaeitootte a ma ea ass 


Pi 
iE 
a 


levet off ‘atalet ©.7 Ate © on): rape ef wm toe 
A Shbeee | 

oh. Gewelt 2h alert wvet ads tveradSGume 1 scgl ig 
oe 


} 

: i " 
_—— » - oe 
f y ay 


20 
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Figure 2.2 A Schematic Diagram of the Apparatus 
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Figure 2.3 Details of the Pressure Chamber 
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2.3 Results 
a) Experiments with Helium: 

A glow discharge was obtained in Helium up to pressures of 300 psi 
gauge (about 22 Atmospheres). The resulting current pulses are shown in 


Figures 2.4 and 2.5 below. 


200 A/cm 
200 ns/cm 


Figure 2.4 Current Pulse Through Helium Discharge at 30 psi Gauge 


The discharge deteriorated at pressures of 350 psi and higher. 
The electrodes used were Aluminium, Rogowski profiled [27] for 1 cm 
separation. The discharge seemed, however, to improve when the gap 


separation was about 2.5 cm. 
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CURRENT 


Figure 2.5 Current Pulse Through Helium Discharge at 300 psi Gauge 


Bigger electrodes (Aluminium, Rogowski profiled for 2.5 cm - 1 inch - 
separation) did not improve the discharge by any noticeable amount. 

The presence of the Strontium 90 sources did not have much effect 
either, and the discharge still deteriorated at pressures in excess of 
300 psi gauge. Six Strontium sources were used in two different positions 
as shown below in Figure 2.6. 

It was noted that when the interelectrode distance was reduced to 


below about 1.5 cm arcing occurred. 
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Figure 2.6 Positions of the Strontium 90 Sources 


Figure 2.7 Tracking Along Chamber sSurtace in Argon 


b) Experiments with Argon 

The above mentioned procedure for Helium was repeated with Argon 
at pressures between 10 torr to 2 Atmospheres. Arcing was observed every 
time. | 

Argon would be very useful in any Hy laser operating as a cid/efi 
laser, because of its high polarizability [14]. 

The arcs, for interelectrode distances larger than 3 cm, (using 
the electrodes designed for a 2.5 cm gap) occurred on the inside of 
the chamber surface facing the 8 sources (See Figure 2d ye 


This was probably due to static charge accumulation because of 8 particle 
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bombardment of the chamber. 


c) Experiments with Hydrogen 

The above experiments were repeated with Hydrogen-Helium mixtures. 
(The minimum partial pressure of Hydrogen experimented with was 5 psi 
while the maximum partial pressure of Helium used was 50 psi). Arcing 


occurred every time. 


d) Experiments with a Ballast Resistor Electrode 
Nineteen 20K2 (5W) ohmite resistors were employed to make the 


resistor electrode shown in Figure 2.8. 
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Figure 2.8 Ballast Resistor Electrode 
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This electrode was used in conjunction with a fine copper wire 
mesh electrode [shown in Figure 2.6(b) above]. The experiment was 


repeated with and without the 8 ray sources. (See Figure 2.9). 


HV 


BALLAST RESISTOR ELECTRODE 


MESH ELECTRODE 


STRONTIUM 90 


Figure 2.9 Position ore Se. Sources 


Without the radioactive sources, it was observed that the presence 


of even a very small amount of Hydrogen (< 1 psi) in a Helium discharge, 


resulted in arcing. 


The radioactive sources behind the mesh cathode did not improve 


the discharge in any significant manner. 
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2.4 Discussion 

The experimental breakdown field strength in Helium is in agree- 
ment (to within 20%) with the value obtained from Raether's Criterion. 
Thus, the breakdown field intensity for He at 300 psi from Figure 2.1 


is about 75KV/cm, while the experimental value obtained is 60 KV/cm 


_ 150KV 
C 225 ean? 


Effect of Dissociation on H, Discharge: 

One reason why it is difficult to obtain a glow discharge in 
Hydrogen is because it shows dissociative attachment [29]. Attachment 
leads to a loss of electrons as carriers since the resulting negative ions 
have a much smaller mobility. It also tends to increase carrier loss 
due to recombination since the cross-section for the latter araceee 
increases. Furthermeres assymetric space charge fields caused by 
attachment can create instabilities leading to a glow-to-are transition. 

Table 2.1 below, shows the dissociation and ionization energies of 
a few diatomic molecules including Hydrogen. It is seen that the ratio 
of dissociation energy to ionization energy is lowest for the Hydrogen 
molecule (28%). Oxygen is the next lowest with a ratio of 0.41 (41%). In 
a discharge, therefore, Hydrogen would probably dissociate first before 
it ionizes. The fact that the Hydrogen atom is very electronegative makes 
matters worse. 

It appears that, although some improvement in the discharge quality 
might have been obtained if a narrower voltage pulse had been used (instead 
of the 200 ns pulses obtainable from the Marx), the major difficulty was 


caused by the effect of dissociative attachment. The strong 8 particle 
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Table 2.1 Ionization and Dissociation Energies for Some Diatomics 


Ionization Dissociation 
Energy, I(eV) Energy, D(eV) 


sources were not capable of making up for the electron losses due to 


this phenomenon. 
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CHAPTER III 


EXPERIMENTS ON RESONANT E-V TRANSFER 


This Chapter describes experiments designed to pump the vibrational 
levels of Ho by resonant transfer from an electronically excited species. 
NO, vapour was excited by radiation from an Argon ion laser. The excited 


molecule was quenched by H, in the experimental chamber. Values for 


2 
the self quenching cross-section and the quenching cross-section with Hy 
are presented. 

It is known that NO, has metastable electronic states (with life- 
times between 55us to 70us) which fluoresce in the red with maximum 
intensity peaks around 6400A and 5900A [34, 35, 36]. These correspond to 
radiation wavenumbers equal to 15,600 a and 16,950 yt respectively. 

The third and fourth vibrational levels of the ground electronic state of 

Hy lie around 14,000 aa and 17,000 ane [18]. It is seen, therefore, 

that the fluorescence from NO, matches closely with these vibrational levels. 
The fact that NO, is a vapour at low pressures (less than 50 torr) and 

room temperature and that it absorbs in the visible [37] makes it a suitable 
material to experiment with. 

Dickens, Linnet and Sovers [38] have considered the quenching of 
an electronically excited atomic species by a diatomic. They considered 
the case where the total interaction Bierce of the final and 
initial states are parallel. Although the quenching cross-sections are 
generally small, a resonance effect is exhibited. Thus the maximum cross- 
section occurs when the energy defect is zero. Dickens et al used a 


simplified model (taking all the electronic states to be spherically 


symmetric). Even though the NO, -Hy system is more complicated the theoretical 
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results of Linnet et al are of general enough validity [39] to justify 


experiments of the type described below. 


(cm) torr-1 x 108) 


A3SORPTION COEF 


WAVELENGTH (A) 


Figure 3.1 Absorption Spectrum of NO, 


3.1 Experimental Investigation of Quenching Cross-Sections 


The radiation from the Ar’ laser consisted of the following six 
° ° ° ° ° ° 
wavelengths: 5145A, 5017A 4965A, 4880A, 4765A and 4579A. Most of the 


laser power (about 200 mW continuous) was concentrated in one 4880A and 
5145A lines. 

The test chamber was evacuated and then filled with NO, at different 
pressures, 

The photomultiplier tube was used in conjunction with a 6180A (red) 


filter to cut off any scattered light from the laser. This ensured that 
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Figure 3.2 Schematic Diagram of the Apparatus 
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Figure 3.3 Enlarged View of Experimental Chamber 


Vyas bh | sy 


~} re 
7 yr met if “IW ' 
Le 
vung 44 
was 7 t. » 
‘ 9 j 


ae oc 


paied =A J de sned alae at oe al 


; ae eae ~ 


*4 


F ool on) ” ti 
7: i a 


Syl 


only a sample of the fluorescence radiation was monitored. 
The output signal from the photomultiplier was measured for various 


and NO,-H, mixtures. 


pressures of NO, g7Hy 


Care was taken to ensure that throughout the experiment the power 
out put from the laser and the pararive positions of the chamber, detector 
and the laser remained constant. This was done by maximizing the power 
output (200 mW) from the laser and adjusting the chopper position (when 
necessary) to get 45 mW average power transmitted past it and by keeping 
the detector flush against the exit window of the test chamber. The 
detector's position was fine adjusted by maximizing the output signal 


from it for a given pressure of NO. (6n6. torr). 


3.2 Results 

Figures 3.4, 3.5 and 3.6 show some typical results for the self 
quenching experiment at low pressures. These show the increasing 
intensity of fluorescence with pressure at low pressures. (See also 
Figure 3.10). 

Figures 3.7, 3.8 and 3.9 show the decrease in intensity with increas- 
ing pressure at higher pressures. For this experiment it was observed 
that the maximum fluorescence intensity, about 1.2 volts, (measured at 
the output of the detector) was constant between about 2 to 7 torr. (See 
Figure 3.10a). 

The amount of scattered radiation picked up by the detector was 
negligible as is seen from Figure 3.12 which was taken with the test chamber 


evacuated. 


Figures 3.11, 3.13, 3.14 and 3.15 show the quenching of NO, meta- 
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stable states by H,- 


3.3 Calculations 


Figure 3.16 shows a simplified energy level diagram of the NO 


2 
molecule with the metastable electronic states lumped together. 
Pumping 
Photon | 
Metastable State | 
| Fluorescent 


Photon 
Ground Level 


Figure 3.16 Simplified Energy Level Diagram for NO, 


Since the on-off time of the chopper was of the order of milliseconds and 
since the longest lived metastable state of NO, has a lifetime of about 
70us, steady state conditions can be assumed in each on-off cycle. 


We may therefore write 
x 
aN*_ go = R - KNAN - = (3.1) 


where 


N* is the density of the metastable state [cm ~] 
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R is the pumping rate of this state fenae sec 


N is the density of the unexcited molecules of NO, iene 


] 

Kis the rate constant for the deactivation collision process fom geen] 

t is the average lifetime of the metastable state [sec.] 
In writing (3.1) the assumption that the other electronic states of NO, 
have lifetimes which are negligible in comparison to the metastable 
lifetime has been made. Furthermore, deactivation by other processes (like 
collision with walls) has been neglected. Equation (3.1) is relevant for 
the self quenching process. 

From (3.1) it follows that 


cu. ee ae Aes 
es (KN +5) 


At low pressures the pumping rate R is approximately proportional to 


the pressure of NO Hence R is proportional to N at low pressures. 


Die 
: , Ne Fivoe ks 
Since the fluorescence intensity is proportional to Fs it is seen from 


(3.2) that at low pressures I, is proportional to N. 


1 


Hence, at low pressures, the fluorescent intensity is proportional 


to pressure. This is the result obtained experimentally as indicated in 


Figure 3.10. 


Self Quenching Cross-Section 
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where C is a proportionality constant dependant on experimental conditions 


(e.g. geometry, detector sensitivity and oscillator strength of the NO, 


metastable states) C would have units of [Joule sec cm]. 


At low pressures the pumping rate R would be approximately proport- 
FO 
ional to the concentration N and hence to 2 , where T is the gas 
KT 
temperature and K the gas constant. Therefore 


2 ee : 
I. i SS (3.4) 
K, ®xo /KT) + (1/t) 1 
2 
Where Ky is another proportionality constant. 
Hence the initial slope of I, vs Px0 (Figure 3.10), m is given by 
2 
CK 
Z —-2 -l 
= Se) 
me oF [J-ciicgebOrr | (3.5) 
Further the "saturation" value of I, is given by 
CK 
1, = a -2 [J om] (3.6) 


; -2 -1 
If S is the sensitivity of the detector in [J cm volt ] then from 


equations (3.5) and (3.6) and Figures (3.10) and (3.11) we get 


CK CK 
ca S and re ics 


Pane Kt 


Hence, K,t = S , from which by using t = 70us [35], T= 300°K, 
19 


= one 
Ko= 1.38 Xx ioe erg/°K = 1.03 x 10 torr cm /°K, we obtain for K, 


the value 
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Converting this into an effective cross-section o, by using [28] 


oT = vo where v is the average molecular velocity at the temperature 


under consideration we obtain 
° 
o = SOE. 


° 
This compares well with the value of 0.9(A)? quoted by Sakurai and 


Broida [40]. 


Cross-Section for Deactivation of NO, by H, 
If quenching by Hy is incorporated into equation (3.1) the following 


equation is obtained for steady state conditions 


NEO ie Ree NEN a Ne 
Toe ee CoD 


where 


K. is the rate constant for deactivation by H,-NO, collisions 


3 BeCaa| and 


[cm 
: -3 
H is the Hydrogen density [cm "]. 


Equation (3.7) may be re-written as follows 
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or 


R 


Pere rnd Kone (3.8) 
1+ Kt + K Nt 


4/2 
+ 


which is the Stern-Volmer relationship [41]. Hence 


I 


ot 1+K,1(P, ,/KT) + Kt (Py /RT) 
Sy paar 5 wy BY RE BY er a a ee Gry) 
vt i+ KU (Pw Kt) K, t (2y/KT) 
where 
Lyre Lav are the fluorescent intensities in mixtures in which the 


Hydrogen partial pressures are Pat and P respectively, 


ja 


Py is the partial pressure of NO, in the mixture. 


Z 


From Figure 3.11 (using intensity values at Pa = 10 torr, Pu = 
13 torr and the value of Ky obtained above) we find the following value 


for K3 


° 
2 
which corresponds to a cross-section value of 0.9 (A) 


3.4 E-V Transfer in H,-NO, Mixtures 


2 2 
To see if vibrational excitation of Hy occurred in collisions with 


electronically excited NO, molecules mixtures containing up to 60 p.s.i. 
(about 4 atmospheres) of Hy aoe about 10 torr of NO, were used. A d.c. 


electric field of about 10 KV was applied to the electrodes. 
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Figure 3.17 Experimental Arrangement 


Vibrational excitation of Hy would have been proved if IR 
radiation were detected. This method of detection of vibrational 
excitation has been used successfully by Polanyi and coworkers [42, 43,44] 


for Hg-CO and Hg-NO systems. 


Since Hy does not have a permanent dipole moment it does not exhibit 
an IR vibrational-rotational spectrum under ordinary conditions. This 
was the reason for using high gas pressures and electric fields in the 


above experiment since under these conditions an induced dipole moment 


could be expected. No IR radiation was observed. 
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3.5. Discussion 

E-V transfer from electronically excited NO, to H, is not an 
efficient process inspite of the near resonant match with two of the 
vibrational levels of Hydrogen. It seems that the quenching of NO, 
observed occurs via an E-T transfer process. 

The results of the above experiment favour the observation that 
there is generally no systematic variation of quenching cross-section 
with the resonant defect [45]. Dickens et aZ [38] have shown that E-V 
transfer, although exhibiting a resonant effect under certain conditions 
(parallel potential curves), would still be highly inefficient for 
processes involving excitation of more than one vibrational quantum. 

This is a result of the extremely small value of the matrix elements 
corresponding to a multiquantum jump. They further conclude that 
efficient transfer may be expected when crossing or near crossing of 
potential curves occurs. 

The present evidence favours the suggestion that efficient E-V 
transfer occurs when a transition activated complex is formed before the 
products. (This will be elaborated upon in the following chapter) It 
appears that a chemical interaction between the reactants, strong 
enough to form a loosely bonded complex in the region of crossing, enhances 
the vibrational energy content of the products considerably [42,43]. 

For the system investigated (NO, -H, ) such a complex, if it were to 


be formed, would necessarily have to be a five atom molecule. It is 


unlikely that such a complex could exist under the conditions in which 


the experiment was conducted. 
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The Landau-Teller model for V-V and V-T transfer [46] in which the 
probability of vibrational transition increases as the ratio of the 
vibrational time period to the interaction time and in which the resonance 
effect is pronounced, is not accurate for the E-V transfer process. 

The activated complex theory requires that the complex have a lifetime 
of at least a few vibrational periods to allow mixing of the electronic 
and vibrational wave functions to occur. 

Nonresonant E-V transfer from the point of view of the statistical 


mechanical theories of Chemical Kinetics is considered in the next chapter. 
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i CHAPTER IV 


NON RESONANT ELECTRONIC TO VIBRATIONAL TRANSFER 


Optical pumping appears to be the most promising method of obtaining 
population inversion in gases like Hy which do not yield to the electric 
discharge or the electron beam methods. This is all the more so since 
to obtain a tunable infrared laser using a homonuclear diatomic (or a 
diatomic having a small dipole moment, like HD) it would be necessary to 
operate the laser at high pressures, as explained previously. 

See E-V transfer has been successfully utilized by Pettersen, 
Wittig and Leone [47, 48, 49] to obtain lasing action in several poly- 
atomic gases. They produced the low lying Br(4°P, 19) state, which has 
an energy of 3685 See by flash photolyzing Br,. This energy level 
matches roughly with the vibrational levels of several polyatomic molecules; 
co, (EOL 021) N,0 (140), HCN (120, 050) and CoH, (01002). This electronic 
state also corresponds to the first vibrational level of Hy (about 4000 
ae but the method cannot be used to obtain an infrared H, laser since 


2 


a chemical reaction occurs in a mixture of H, and Bry yielding HBr [50, 
Slee oat. 

Fortunately, a resonant match between the levels is not an essential 
prerequisite in many other efficient E-V transfer processes. Such 
processes have long been known to occur [52.53.54] Here, .elticitent 
conversion of electronic to vibrational energy occurs essentially because 
of a mixing of the wavefunctions representing these types of energies. 


A loosely bonded activated complex molecule is formed which then dissociates 


into vibrationally excited products. 


48 


gM hs irae: 


qrininas he Nokpda sndebauny Peapankay 7 Dee Oe 
cool to ie lh al zs s mrNeY 
- eueie oa psn a4 Iipses chi. y 
e 26) sintected’ accion a sina aid ale io 
on Tagen ot bingw thy iodl ott aan 
own sven Reehalorce 4a. yi . 
oreo iet xd beet lity et cedosie ned el 
-qieg f.powte mw cotpen geet Mantowe (oe es 
col Uptiw were ri al guys wok wp nese ea oegs 
feet ae wei pe anaprentale: dwell 2 @ . - a 
whiole ordaetay Laven So: eter! (ynmdtetdey: 1 anv gages | 
ainevinec Le atat Di id ray ee be (ome UA i mes | Ges 


Gas tredal Ji 30 \Gaeed dene sandy qq aas the 
ewts sexy gi Sota Rah ae manedn we haw ee sans whee ts neh 
OEY eke gab Tay. bla * ali aia. paaa a 
ea 7 Lied el 

atoms at an sl 9 ent aR, 
Aye tepony Sain: vay ‘tahiti a ta a ape vod vi 
naniriht yo [OR fhe St) wwsdeee wernt vend mr wed pads is 
went yhhst tiene cwaae, iectain, (apebiatty oe" Salinas n oxieaa hit 
witysis ~ 4s jeans) gis) 2 -isgrinsei toll Ww. patie a w- 4 
> abopmes were ae heme wh ain enki Sealife | aie: coal Chet yas nonl he isa 


a > ceatnbeeny bio od ie fp ty ay : iw LL? crak a fj 


| 


F 1 
ee. * al , : 2 : Toe and 


49 


The Activated Complex Theory of chemical kinetics, and its extension 
by Rice, Cen ne cer: Kassel and Marcus (the RRKM Theory) are general 
theories which apply to many chemical processes in which energy re- 
distribution occurs [56-59]. Some processes like the quenching of elect- 
ronically excited oxygen (o!D) atoms by N, in which no new chemical products 
are formed can also be treated by this theory [60]. 

For reactions on which enough data are available, the RRKM Theory 
allows a quantitative estimation of the amounts of vibrational and trans- 
lational energy content in the products. zor instance, Tully [60] has, 
by an extension of the RRKM Theory, successfully explained the large 
cross-section in the quenching of o(tp) by N 


2? 


+ O(°D) +N, + AE 


op) +N ; 


2 

This Chapter is divided into three main parts. Section 1 examines 
the activated complex and the RRKM Theories. Some relevant background 
material is also presented. In Section 2, a calculation for the cross- 
section of quenching in the Na (7P) + Hp > Na (7s) + H, + AE reaction is 
given, together with the (partial)cross-section for the excitation of the 
second vibrational level of Hy « 

Finally, the E-V transfer process is viewed in terms of the break- 
down of the Born-Oppenheimer Approximation. It is noted that under these 
conditions, since the electronic and nuclear (vibrational and rotational) 
motions are no longer separable, efficient E-V transfer may be expected. 
From this point of view, the necessity of a loosely bonded activated 


complex [54] and non-parallel potential curves [39, 63] follow automatically. 
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“he conditions under which the Born-Oppenheimer Approximation breaks 

down can also be used as a basis for selecting atoms (or molecules) which 
would pump (by E-V transfer) a potential laser molecule. A general 
criterion for this choice is meeaenteds 

ae Chemical Kinetics [56, 57. 58] 

It must be noted that there exists a great nent of confusion in 
the current literature regarding notation in this field. Robinson and 
Holbrook [57] have attempted to rectify the situation by defining a 
consistent set of symbols and comparing them with those of other workers. 
Unless otherwise stated, their notation will be followed here (see 
Appendix V). | 


Consider a reaction of the following type 
D-+(B -— C)>(D — B+ C (4.1) 


in which a diatomic molecule B-C reacts with an atom D yielding D-B 
and the atom C. 

Many such reactions can be analyzed in terms of a potential energy 
diagram like the one shown in Figure 4.1. 

The reaction path shown in Figure 4.1 is the path of minimum energy 
and is hence the most probable one followed. The "saddle point" or 
the col shown on this path has the highest energy on the reaction path. 
Any other path would have a maximum at a higher energy level than the 


saddle point. 


A section through the reaction path is shown in Figure 4.2. 
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Figure 4.1 A "typical" potential energy surface diagram 
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Figure 4.2 Section through a reaction path 
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The Activated Complex Theory also known as the Absolute Rate Theory 
and the Transition State Theory, postulates that a loosely bonded molecule 
exists around the region of the highest point in the reaction path.© this 
molecule is called the activated complex and is in equilibrium with the 
"reactants', D and B-C. 


Hence, equation (4.1) is more accurately written as 
D +(B - C)z D°--B - C +(D — B)+C (4.2) 


The theory stipulates that all the activated complexes dissociate 
into products. Hence, although the probability of an activated complex 
being formed is, in general not wate once it is formed the Activated 
Complex Theory assumes that the probability of passage over the potential 
barrier is unity. 

Recent work (theoretical and experimental) has shown that, in general, 
the reacting particles form an "Energised" Molecule (collectively denoted 
by A’) which has enough energy to dissociate into the products. However, 
Me is distinct from the Activated Complex in that although the former has 
enough energy, the distribution of this energy may not be conducive to an 
immediate dissociation into products. It is only after this molecule 
takes on the right configuration, in which the energy is favourably dist- 
ributed in its different degrees of freedom, that the reaction goes to 
completion. For this reason, some authors prefer to call the Activated 


Complex by the term "Critical Configuration" [59, 60]. (In our notation 


os 
the Activated Complex will be denoted by A ). 
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Equation (4.2) represents a chemical reaction. However, the 
Activated Complex Theory applies to many processes, in which only the 
energy in participating particles is redistributed and which cannot strictly 
be called chemical reactions. An inane eee of such processes that 
can, nonetheless be treated by the Activated Complex Theory and its more 
accurate extension - the RRKM Theory - is the quenching of electronically 
excited atoms by diatomic or polyatomic molecules. 

The existence of activated complexes has been proven experimentally 
in many instances [61]. Furthermore, the concept of a reaction path 
easily explains why for certain reactions the products take up most of 


the energy of exothermicity as translational energy, whereas, in others, 


the conversion is essentially into vibration. This is illustrated in 


Figures 4.3 and 4.4. 
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Figure 4.3 A potential energy surface of the "attractive" 


type, for the reaction D+BC+DB+C [56] 
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Figure 4.4 A potential energy surtace ot the "repulsive" 


type, for the reaction D+BC+DB+C [56] 


Figure 4.3 illustrates the type of reaction where the molecule DB 
might be expected to take up most of the energy released in the form of 
vibration since the activated compen occurs at a position where the D-B 
bond is stretched. Such a potential energy surface is sometimes called 
the "attractive type" as opposed to the "repulsive type" shown in Figure 
4.4. In the latter, the activated complex occurs at relatively large 
distances of C from DB. Hence, most of the energy released goes into 
translational kinetic energy of the products. In general, the potential 
energy curves are such that the products absorb the energy of reaction in 
both these modes although the distribution is usually uneven. 


Quenching reactions are quite similar although no new chemical products 
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occur. | Quenching may be expected to produce either vibrational or 
translational energy depending on the energy distribution at the critical 
configuration. 

The potential energy surface concept holds, in principle, for any 
activated Bancless An N-atom complex molecule has 3N-6 vibrational degrees 
of freedom if it is nonlinear and 3N-5 vibrational degrees of freedom 
if linear [62]. Hence, a complete potential energy plot for such a molecule 
would be an n dimensional plot (where n = 3N - 6 or 3N- 5). This being 
hard to draw, the procedure is to hold some of the coordinates fixed and 
depict the variation of potential energy with the rest of the coordinates. 
Hence, Figures 4.3 and 4.4 are valid plots for a nonlinear triatomic system 
with the angle between its two bonds fixed. It could also represent the 
plot for a linear triatomic with the two bending vibrations' normal 
coordinates held fixed [57]. 

(It might also be pointed out that whenever a potential energy 
diagram for a molecule is plotted,it is implied that the vibrational motion 
of the molecule takes place under the influence of this potential, for 
a given electronic state [18]. This, in turn, implies that the nuclear and 
electronic motions of the molecule can be separated by the Born-Oppenheimer 
Approximation. Under the condition when this approximation is not accurate, 
the potential energy surface concept, strictly speaking, does not hold [18]). 

Consider the process indicated by (4.2). The activated complex theory 
makes the hypothesis that the activated complex (D -+-B-C) which we will 
denote by Ne is in equilibrium with the reactants D and B-C even when the 
reactants and the products are not in equilibrium amongst themselves [56]. 


Under these conditions, statistical mechanical theory [56, 57] leads to 
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the following expression for the cverall rate constant k 


d 
ee + 
dt KTe-Q : 
k = SO = exp(-E_/KT) (453) 
[D][BC] “h QQ, 0 
+ 
where Q , Qp> Qa are the partition functions of the activated complex, 


the atom D and the molecule BC (in evaluating Om the activated complex 

is considered to have one less degree of freedom than would be normally 
expected, since the derivation of equation (4.3) treats the motion of 

the complex over the potential barrier as a degree separate from the others; 
the effect of this is expressed by the factor = in the above equation). 

E. is the difference between the zero point energies of the activated 
a ie and the reactants, which may be regarded as the activation energy 
for the reaction at the absolute zero of temperature (see also Figure 4.5). 

k and h are the Boltzman and Planck constants respectively, and 

T is the absolute temperature (See Appendices V and VI) 

The Activated Complex Theory has been successful in explaining reaction 
rates for see processes. At low pressures however, there is a substantial 
discrepancy between the experimental rates observed and the theoretical 
ones predicted. 

The RRKM Theory is an extension of the activated complex which has 
found the most success from amongst Pi oot day theories. It postulates 
the existence of an energized molecule before the activated complex is formed. 
The Activated Complex Theory, on the other hand, does not distinguish 
between these two states. 

Before going into the mathematical formulation of the RRKM Theory, 


it might be useful to clarify some of the terminology commonly employed. 


For this we refer to Figure 4.5 and to Appendix 5. 
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Figure 4.5 Energy level diagram relevant to the RRKM treatment [56,57] 


An energy mode is said to be adtabatte if it stays in the same 
quantum state during the reaction. Generally, an adiabatic mode has 
the same energy content throughout the reaction. novevee. this may not 
be true for an adiabatic rotational mode since in many instances the 
configuration of the system and hence its moment of inertia changes 
considerably during a reaction. Hence, although the rotational state may 
be characterized by the same quantum numbers, its energy content might 


change. 


Adiabatic rotations are denoted by Ey and Ey in Figure 4.5. 

Any form of energy which cannot be redistributed is called fixed 
energy. Examples of fixed energies are the zero point levels of molecular 
vibrations, E. and Eo in Figure 4.5. 


The RRKM Theory as normally applied assumes that all nonadiabatic 


degrees of freedom can take part in the redistribution of energy in the 
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various modes of the system and are hence called also active modes. 

An energized lee ae. is defined as a molecule which contains, 
in its active degrees of freedom, a nontined energy En greater than Eo 
(the difference between the zero point levels of the reactant's and the 
activated complex). ES is the activation energy of the reaction at absolute 
zero of temperature and is hence, classically, the critical value of a 
below which reaction cannot Baers 

The activated complex is a species having a configuration corresponding 
to a position at the top of the energy barrier in the reaction path. 

The total active energy in the energized molecule denoted by a can 


Send k 
be resolved into a vibrational component, EY and a rotary component, EN. 


where 


ae oe and E. are the vibrational, rotational and translational 
Veer 


+ 
(along the reaction coordinate) components of E . 


The RRKM reaction scheme is 
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Equation (4.4a) denotes the formation and deactivation of energized 
* 
molecules A ., 
Equations (4.4b) and (4.4c) show the formation and dissociation of 
the complexes Cys respectively. 
k)> ko and kr are functions of the corresponding energies, this being 
another major difference between the RRKM and Activated Complex Theories. 
A , : A 
OK, (BE is SE*) is the rate constant for the reaction which produces 
co * * * 
energized molecules having energies between E and E + 6E. 
Applying the steady state condition to reactions (4.4a) and (4.4b) 


we obtain the following equation, in which square brackets denote the 


concentration of the appropriate species, 


* * * 
= 6k [A] [M] - k,[A ][M] - k,(E )[A ] (4.5) 


* ok * ; 
(The subscripts (E> E + SE ) have been dropped from equation (4.5) for 


convenience. Hence, 
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The overall rate of reaction (formation of activated complexes - 


which is the rate at which particles arrive at the critical configuration) 
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cee ; 
Ane ke ) 6k, /k, [A] 
k,(E [A ] = ~——_,-—+—_ (4.7) 
1+ k € )/k, [M] 


The unimolecular rate constant for the process is defined by the 


following equation 


Peet 
Kuni ~ [Ay ( dt ) 
or 
afA] _ 3 
Sages Toni | (4.8) 


From equations (4.7) and (4.8) the overall unimolecular rate constant 
for the processes denoted by equations (4.4a) to (4.4c) is given by 
k (E')/k 
a ) Z 


Ser oe SET ey 
1+ Kk, (E")/k, [Ml] (EE der) (4.9) 


fe) 


The lower limit in the integral is Eo because this is the difference between 
the zero point levels of the reactants and activated complex molecules and 
hence is the minimum energy required for the formation of the latter. 

To evaluate (4.9), we need an expression for k(E), the rate constant 
for the formation of activated complexes from energized molecules having 

* 

active energy E. This expression is also fundamental to any analysis 
of the vibrational energy of the products in which we are interested. (This 
will be seen below in Section 2). 


In the steady state, we find from equations (4.4b) and (4.4c) 
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(a) CE) = CER (4.10) 
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> + 
where A denotes the activated complexes moving in the forward direction 
(towards dissociation into products)only, over the saddle point in the 


reaction path. In most cases [57] it can be shown that to a good approx- 


imation 


Sey ; : : ; ; 
where [A ] is the total concentration (including complexes moving in 
both directions). 
Hence 


- 
er elie ae tA 
ke ) = 5 k CD area Gilt) 


[A ] 
Since the active energy of a complex can be divided, in different 
d ; P an 
ways, into the energy of vibration and rotation (EO. and the energy 
: ‘ + 
of translation in the reaction coordinate (E,)> the processes denoted by 


(4.4b) and (4.4c) can, in effect, be resolved into those represented by 


(4.4d) below. 


+ + + + 

A te a K(E,) 

A (EF ) ———>> A gt zt) ——— > Products (4.4d) 
aes is 


where we have used the fact that 


* 
E = ET om 
fe) 
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The rate constant k (E ) can hence be written in the more general form 
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+ ete 

KE ana) A Gaeta ae (4.12) 
at A [A (E*) ] 
NAG 


To find K (EL), the decomposition of the complex into products may 
be treated as the cranaiatien of a particle of mass ‘ in a one-dimensional 
box of length 5. where the latter is an eehdrrary (small) region at the 
top of the energy barrier which may be Bail to define the activated complex. 
(We note that the final expression for k, (E") is independent of both 


u and 6. See below). The speed of translation is (Et su)? /? 


{2 


antes ces 
time taken to cross the barrier is hence 6/(2 E./u) - The rate constant 


and the 


is the inverse of this quantity 


K’ (E*) = (Et /us?)1/? (4.13) 


* 
The ratio [AT J/[A ] has to be evaluated for the derivation 
+ + x 
(ESE) (E ) 


Wag AE 
k 
of kG ). It should be noted that in equation (4.12) this ratio is for 
Se % bes . 
the values of A and A at equiltbriun. 


Statistical mechanical considerations [57] enable us to write this 


ratio as a ratio of the corresponding partition functions. 


Hence, 
Me fi 
+ + ee ie 
(t.) neg (at yin (Deg ke EL) | anes 
Acgty Q(a’) (Je, exp(-E /KT) 


eqn 


where 


t ‘ £ 4 #8) ‘or 
i 


1-7 ikea areal & 


ryA Ge a Tr yaw 60°. Grae ou ‘eh ) ea? et ; 
‘ al CT het ‘ ad ‘ So 3% sad Lith iyo ci vo. ; pepe af 
hes hei Szyk, alte dah er Tadetiaite al Pigg ta iat 


er wae ihe ch Mier at aaa ebhneite ee . 
| ’ hunt ae wo 


. 


j ; iyvsiy 6679 Wha. saad ou Sti | f 5 A at t iy m | 2 
| f 7 Bat-\k ob ats At ity Paks ' 
onsen ety & ‘een ‘oa inv mal 7a 


views Os bite ana) oe ra lavays Boge Aton 


YAR 


63 


Q(A*) and Q(A*) are the partition Eugerions of the activated 
complex and the enerei el molecules respectively, 

leq and Je, are the number of Geant pratee in the energy range 
E46 = E+ 6E 


Equation (4.14) therefore reduces to 


e bb... 
oe NC be) | 
—— Siler ae (4515) 
[A ] Ney 
eqn 
where the subscripts have been dropped for convenience. 

ae ee Epa ialchey ob + + + 

NE ey Be) ae oe voce (4.16) 


where 
te : A ' : ate 
P(E. is the number of vibrational-rotational levels of A 
i ae 
with vibrational-rotational nonfixed energy equal to Er 
and 


Nv (E,)5E, is the number of states with translational energy in 


+ 


+ 
! os SE, 1571s 


+ 
the range E. aE 
Quantum mechanical treatment [65] of the translation of a particle, 
mass u, in a box of length 6 produces the result that the energy is 


z ar 
quantized and n, the number of levels having energy up to and including E. 


is given by 


(4517) 
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h is Planck's constant. 


The density Nv (E,) is therefore given by 


Y/2 


dn a ( ous? ) 
+ 2_+ (4.18) 
dE h Ee 


gee ce 
ee _ 


Substituting (4.13), (4.16) and (4.18) into (4.12) the following 


important RRKM expression is finally obtained 


hee 
kK )=— +, J pat) (4.19) 
aes te 
E =o 
vr 


As might be expected, when equation (4.19) is substituted-into (4.9), 
the RRKM expression for K nd reduces to the Activated Complex Theory 
expression (4.3) in the high pressure limit ([M] > ~) [57]. 

Equation (4.19) was first derived by Marcus [59] and is the starting 
point in the calculation (section 2 below) for the vibrational energy 
content in a quenching molecule after a collision with an excited electronic 


species. 


4,2 E-V Transfer 


Many quenching processes of the type 
* 
D +AB> D + AB(v) + AE (4.1a) 
where 


* 
D is an electronically excited atom, 


AB(v) is a vibrationally excited molecule, 
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and 
AE represents the energy defect, can be investigated by extending 
the RRKM treatment described above. 
Tully [60] has applied the fundamental RRKM expression, equation 
(4.19), to the quenching of electronically excited Oxygen atoms by 
Nitrogen molecules: 


OC Dene 0 ry N,(v) + AE (4.20) 


7 
The treatment developed by Tully is especially for spin forbidden reactions 
like (4.20). The principle, however, is peneret enough to be applied to 
other reactions too. 

The RRKM theory has been particularly successful in explaining the 
high efficiency of the E-V transfer process in reaction (4.20) which 
cannot be explained otherwise. Thus, for instance, the Wigner Spin 
Conservation Rule predicts that a reaction like (4.20) should occur only 
with an extremely small probability (being "forbidden" under the Rule). 
However, it has been found experimentally that the reaction has one of 
the highest rates observed for any process of this type. The RRKM Theory 
explains this otherwise anomalous result by postulating the formation of 
a relatively long lived energized molecule (0°). 

Equation (4.19) is rewritten to include the energy of the activated 
complex and the angular momentum explicitly. kK is also assumed to 


be dependent on the translational energy of the complex. The equation 
therefore becomes 
+ 
aoe 8) =e) (4.19a) 
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where 

mut ake : 

E is the total active energy of the complex 

+ + + +, t,t, ot 

= B, + +E =B, +E + ; B, bei iti 

(E B, Er Ee B. EY ED + ELS B. being an additional 
form of active energy defined here as the height of the centrifugal 
barrier that the activated complex has to cross before dissociating 
into products. Note that in the usual RRKM treatment, ET is 


+ a 
written as E =E 


~ ~ ~ + 
+E =E +E _+E,. Furthermore, the sub- 7 
vr & Vv r t 


script i is used to emphasise the fact that there may be more than 
one reaction path leading to the formation of an activated complex), 
J is the total angular momentum of the complex and 


+ + 
WCE.) denotes the sum }) P(E.) 


Another improvement is made by assuming that all activated complexes 
do not necessarily decay into products as the normal RRKM theory assumes. 


+ 
E 


ae B,), is included on the right 


+ 
Hence, a probability factor T(E,» 
hand side of equation (4.19a). 
Equation (4.19a)is therefore modified to equation (4.19b) 
=e (4.19b) 
hN (E ) 


4.2.1 Quenching Cross-Sections 


Equation (4.19a)may be viewed as describing the rate of arrival 
at the critical configuration, while (4.19b) describes the rate at which 
these activated complexes emerge from the critical configuration leading 
to products with translational energy E, along the reaction coordinate. 
Tl. is assumed to be close to unity for reaction paths not involving 


curve crossing. For paths that do involve crossing (or "avoided" crossing) 
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Ty is assumed to be given by the Landau-Zener [66-69] formula for 
transmission coefficient: 


V7 


; im shen re ats 
E> B,) = ,(B) = 1 - exp{-a,/(B) '7) (4.21) 


A, is the Landau-Zener coupling coefficient for the electronic 
transition under consideration. 
Note that in equation (4.21) we have made the assumption that the 
ve = Sr A 
dependence of Tr, on EY is stronger than that on Er and Bi. This 
seems reasonable if it is borne in mind that Pr. reflects the probability 
of an electronic transition occurring at critical configurations that are 


defined by curve crossing. 


Activated 
Complex 1 


Figure 4.6 Effective potential along some arbitrary reaction coordinate 


illustrating the two paths by which critical configurations may be reached [60]. 
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+ 
Tf A (E »J) is defined as 
+ 
E -B; (J) 
+ + He oe 
A(E ,J3) = f RACE 4 s35E. YdE (4.22) 
ai eG at ie t 
then it follows that the probability that an activated complex, nee with 
ee 
total active energy E and angular momentum J, will decay into products 
via channel i is given by 


+ 
A, (E J) 


6,(E',J) = (4.23) 


yi he Cees) 
j J) 


where the summation is over all possible reaction paths leading to an. 
activated complex. 
The centrifugal barrier height B, (J) is needed in order to 
evaluate the definite integral in equation (4.22). 
This can be found if the form of the potential energy (in Figure 
4.6) around the critical configuration is known. Let this value of potential 
energy be denoted by Me | 
Around the critical configuration on the reaction path the total 
energy (or "the effective potential barrier") [18] will be given by 
Te, Olay Ge) oy = (4.24) 
eff if 
where 


u, is the reduced mass of the complex 


+h 


R, is its effective radius 


he second term on the right hand side of (4.24) is just the 


Hoe 
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kinetic energy of rotation of the complex, hence the name centrifugal 
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barrier for V oe¢) = B,(J)]. 
max 
It will be assumed that for channels not involving curve crossing 


Phot valuelor v; witli be -Gemthe: form 
ob om 
V,(r) =-Cr 6 (4525) 


where 
C is the van der Waal's coupling coefficient for the interacting 
particles, 
r is the effective distance between these particles. 


Using (4.25) in (4.24) differentiating with respect to r and setting 


ay 0 to find the maximum yields 
cy 
ES eae 
J 


from which 
4 eet ae 


This is the value of the centrifugal potential barrier for a 
critical configuration not involving curve crossing. 
For the critical configuration(s) involving curve crossing, equation 


(4.24) is used. This requires that the coordinates of the points of 


intersection be known reasonably accurately. (As will be seen below, this 


restriction usually causes difficulties since potential energy surfaces 


for many processes are known only very roughly, if at all). 


The cross-section for the reaction occurring (to completion) via 
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channel i is given by 


max a 
2, ORL ‘les e(J) $,(E ,J)ds (4.27) 


where 
E. is the initial collision energy, which is a function of the 
translational temperature, 
AL is the cross-section for complex formation, 
o(J) is the normalized distribution of angular momentum of the complexes. 


eee can be written in the form 


2 
Ok Cl aes (E,)/2u,E, (4.28) 


where 

Bs is the statistical degeneracy of channel i. 

This follows from the fact that the effective collision cross-section 
is directly proportional to the kinetic energy due to rotation and in- 
versely proportional to the centrifugal barrier height. 

To evaluate 0, (E) in equation (4.27), the values of the total 


4. 
E 
. + ~ 
number of active vibrational-rotational levels, WE) = ) P(E in 


the activated complex, and the number of energy levels of the energized 
molecule per unit energy, NCE), have to be obtained. 

This is done by treating the complex as a set of S off classical 
harmonic oscillators, each oscillating at a frequency V5 [57]. | 


The value of Soft is usually taken to be 
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+ 
See = (4.29) 


where 


s is the number of vibrational modes at the critical configuration, 
r is the number of active rotations of the complex [57] 


The following expressions result 


A Pate 
(E" ) ants 
W(E te ee (4.30) 
1 
vr (Scge? I (hv, ) 
-1 
and ‘ Soff 
ko (E ) 
Ne aa ae (4.31) 
(Sigg Ty! (hv, ) 
where S 
eff 
I(hv,) is the product II (iva. 
ak j=1 SL 


(The derivations of these expressions are given by Robinson and Holbrook 
in reference [57]). 

If a system having only two paths leading to critical configurations 
is considered (path no. 2, involving an intersection of potential curves 
and the other, path no. 1, not) then we obtain the following expressions 


[57] by substituting equations (4.30), (4.31), (4.21) and (4.19b) into 


equation (4.27) 


en max A) 
a) oun, ; ru vanes (4.32) 
7 Snax A’, 
errr anrere (4.33) 
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where 
0, (E,); 0, (E,) are the reaction cross-sections for paths 1 and 2 
respectively, 
at, = = [EY = 27 (54cu3y/? 7? 
and 
A’, = 24,[E - TGicn e121 A, being the Landau—Zener 


coefficient for path 2. 


In deriving (4.32) and (4.33) the value of S, has been taken to 


££ 
be 2. This is the "loose complex" limit for a triatomic activated complex 


[60]. Furthermore, the normalized angular momentum distribution has been 


taken to be [60] 


2 
0 (J) 2I/I Forma) %< Lie (4.34a) 


x x 


and 


i] 
j=) 


o(J) for J > sey (4.34b) 


x 


The value of a ae is the value of angular momentum corresponding 
to the maximum value of V of (the centrifugal potential). From (4.26) 


it follows that 


_ 2 SLi 4.35 
at Me (S4E Cu) (4.35) 


The total cross-section On (E,) for the reaction (denoted by 


equation (4.la) is hence 
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Op (E,) = 9, (E,) + 9, (E,) (4.36) 
4.2.2 Vibrational Energy Content of Products 

An estimate of the product vibrational energy is essential for the 
analysis of any infrared diatomic laser, since such a laser would require 
a partial vibrational-rotational population inversion. It is therefore 
necessary to calculate the extent to which the exothermicity of the 
quenching process is partitioned into product vibration. 

‘ In order to estimate the extent of vibrational excitation in the 
products formed by dissociation of the activated complex, it is necessary 
to make certain assumptions regarding how the energy of the complex 
evolves into translational and internal motion of the products. 

Tully [60] makes the assumption that the total rotational energy 
of products equals rotational energy at critical configuration, vibrational 
energy of products equals vibrational energy at critical configuration 
and translational energy of products equals the sum of the translational 
energy at critical configuration plus the effective potential energy EBORE 
This leads [60] to the following expression for the partial cross- 
section for production of vibrational energy in the diatomic molecule (AB 


in (4.1la) above) 


+ + + 
Bch.) Es 


mex 25 1 f 1/2 


[1-exp(-A,/E,)"“)] 


dE, aJ (4.37) 


where 
E is the vibrational energy of the diatomic. 


(The other symbols have been defined previously). 


Jessinay., s, & 


ae ; ine 4 ha on 
it = he 


tals Oh Daeg Ha od ti + eid ssc 


18) aoe Oe? oe? Ep eee a ee 


gibys. i. Prey eit ay sara 


aks oleae Foy ehier el any nh poe 


as 


a 
oe 7 


eer a a 


pris io Soi a ‘opie any at va 
: ee a a ated ial 10st op 


74 


k 
Application to the Na “H, quenching process 


A computer program utilizing the Romberg technique of integration 
[70] has been developed to evaluate integrals of the types shown in equations 
(4.32), (4.33) and (4.37) above. (Appendix VII gives details of the 
program). This routine was tested by duplicating the results of Tully [60] 
for quenching of electronically excited Oxygen atoms by Nitrogen molecules 


OD) eee 210 CE) 4H N, (v) (4.38) 


2 
(The theoretical results obtained by using Tully's analysis predict 
an efficiency of about 30% for this Reyeceee ey process, in very good 
agreement with the experimental results of Slanger and Black [71] who 
report that an average of 33% of the 1.98 ev exothermicity goes into N, 
vibrations). 
The reaction denoted by equation (4.39) below which results in 


vibrational excitation of Hydrogen molecules has been known for quite some 


time [52,39] 


> Na (378) + H, Vey (4.39) 
Vv v* 


Na (3°P) - Hp 


where 


v and v' are the initial and final vibrational quantum numbers. 
Quenching of the excited Sodium atoms occurs via the formation of 
a loosely bonded ionic complex, NaH, the only restriction being that 
the Hydrogen bond be extended initially [52].. (Hence, in (4.39), v # 0). 


The RRKM theory of unimolecular decomposition described above can, 


therefore be applied to this process. 
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The reaction indicated by equation (4.39) is not spin-forbidden, 
in contrast to the quenching of excited Oxygen atoms by Nitrogen molecules 
(4.38). Therefore, Tully's extension of the RRKM theory would be expected 
to yield only a rough estimate for the vibrational energy content of H,- 
The use of the theory can, however, be justified on the grounds that 
Tully's formulation is fairly insensitive to the Landau-Zener coupling 
coefficient A, (this being the reason why the approximate Landau-Zener 
expression could be used in the first place) and because the most important 
criteria are the formation of an activated complex and an intersection of 
potential Bree, curves. Both these criteria are satisfied by reaction 
(4539): | 


The relevant potential energy curves are shown in Figure 4.7. 


Oo 


0 Na H_ distance 


Figure 4.7 Potential Energy Curves for the Na-H, System (assuming linear 


Configuration [52] 
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Magee and Ri [72] have shown that the activated complex in process 


(4.39) is actually triangular in configuration (see Figure 4.8) 


20; 


Figure 4.8 Configuration of the Na"H Complex [72] 


The potential energy curves in Figure 4.7 are, strictly speaking, 
valid for a linear configuration. However, the error induced as a result 
of this assumption is negligible [52]. 

To be able to evaluate the total quenching cross-section (equation 
(4.27)) and the partial vibrational cross-section (equation (4.37)), 


the following datawere used: 
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(See Appendix VIII for a note on units). 


“ay Na 
WPStgoea yee is the reduced mass of the activated complex, 
H, Na 
where M. , M. denote the masses of H, and Na,respectively. Using M., = 2 a.m.u 
H, Na 2 Hy 
M. = 23 a.m.u. yields p = 1.84 a.m.u. 


Na 


C is the van der Waals' coefficient for the Na-H, system. This can be 
calculated from the combination rule for van der Waals' force constants [73]. 


This rule gives for a system a-b of particles the following value: 


2 CAs Sat 


oer (a, /a dC. (a /o Cy (4.40) 


where 
the C's are the induced-dipole-induced-dipole van der Waals 


coefficients for pairwise interactions of atoms a and b, 


and 
the a's are the static dipole polarizabilities. 
Using the values given by Kramer and Herschbach [73] equation (4.40) 
yields 
C ps So Sls dle x00. 90"x Ot erg nw 
Na-H, 3:33, 1580 + 222. x 13.0 
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Z 4 
hen ~ 51.2 ev(A) 


EY is the exothermicity of the reaction 
+ 
E = 2.09 ev [39] 
+ a ; : ‘ ‘ 
The values of Rand V, are obtained from Figure 4.7 which is taken 
from reference [52]. 


The value of the Landau-Zener coupling coefficient, A is uncertain 


2? 
for the quenching of Sodium by Hydrogen. However, as pointed out previously, 
the formulation given above, is not very sensitive to variations in Ay. 
The value of Ay used here is taken from Yamanouchi and Horie [74] and is 
the value for quenching of excited Oxygen by Nitrogen molecules. 

The graph in Figure 4.9 summarizes the results of the computations 
done. It is seen that at the lower collisional energies, Ew the per- 
centage of exothermicity which goes into vibration is larger than that at 
higher temperatures. The total quenching cross-sections show a similar 
trend. This behaviour is to be expected because at very high temperatures 
the duration of the collision becomes so short that the particles do not 
interact by any significant amount. Note that the calculated values of 
vibrational cross-section [= 10(A)7] for the excitation of the second 
vibrational level are in fair agreement with the experimental value 
[= 15(A)7] quoted by Pringsheim [79]. [The latter value is for excitation 
of the fourth vibrational level. However, since the RRKM theory for E-V 
transfer predicts that vibrational cross-sections are not heavily dependent 


on the resonant defect - and this has been verified experimentally [39] - 


we would expect the experimental cross-section for excitation of the second 


ony: 
level to be about 15(A) too]. 
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Figure 4.9 Computed total and partial-vibrational cross-sections for 
quenching of Na? by H,- 
4.3 The Born-Oppenheimer Approximation and Electronic-to-Vibrational 
Energy Transfer 
The Born-Oppenheimer Approximation allows the total wavefunction 
of a molecular system to be separated into the product of nuclear and 


electronic wavefunctions. The total energy of the molecule is also 


79 


separable, by this approximation, into the sum of nuclear and electronic 


energies. 
The importance of this approximation to the E-V transfer problem 


lies in the fact that if, in an interacting system of particles, the 


nuclear and electronic motions are not separable then it might be expected 


that the E-V and V-E processes would occur with high efficiency. The 
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breakdown of the approximation and the consequent mixing of Bieeeronte and 
nuclear wavefunctions allows efficient conversion of energy from electronic 
to nuclear (vibrational and rotational) forms and vice versa. 

This statement is verified by the fact that interconversion of 
electronic and vibrational energies occurs with high efficiencies in systems 
where the potential energy curves of the final and initial states cross 
with formation of a collision complex at the intersection. This is 
exactly the condition where the Born-Oppenheimer Resear ime ction breaks 
down, as will be shown below. 

This Approximation was first derived in 1927 [75]. The original 
derivation is fairly complex [76] and there has been little work done since 
then in this field. However, the validity of the approximation has been 
proved theoretically, for a simple system recently [77]. The Approximation 
has been verified, at least indirectly, by experimental evidence. Hence, 
for example, since the Franck-Condon Principle [76] is based on this 
approximation and since many molecular phenomena have been explained by 
the former, it may be said that the Born-Oppenheimer Approximation has 
been proved experimentally. The E-V transfer frocess provides indirect 
experimental evidence for the validity of the approximation, also as explained 
above. 

It is, perhaps, interesting to note that the Franck-Condon Principle 
rests on the validity of the Born-Oppenheimer Approximation. Efficient 


E-V transfer, on the other hand, relies on the breakdown of the approximation. 


Derivation of the Approximation [78] 


Classically, the picture that forms the conceptual basis of the 
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Born-Oppenheimer Approximation is one in which electrons necae tenatieed 
as light and extremely swift particles "whizzing furiously" around heavy, 
and sluggish nuclei. It might be expected that the electrons complete 
traversing one or more orbits before the nuclei move even a very small 
distance. This model is based on the fact that nuclei are immensely more 
massive than electrons. The same Coulomb force would consequently give 
an electron a much higher acceleration (and hence a much higher terminal 
speed) than it would a nucleus. 

To derive the approximation quantitatively, however, we have to 
resort to quantum mechanics. 


The starting point is the Schroedinger time independent equation 


where 


a> Ry denote the radius vectors to the electron i, and the 


nucleus, 4a, 


E denotes the energy of the stationary state described by equation (4.41), 


H is the Hamiltonian Operator for the system. 


2 e a 2 
=e ise h p34 h 2 tat Fes 
BES hn re Ue are ye Ear (4-42) 
fe) a 
The first two terms in equation (4.42) are the Kinetic Energy operators 
for the N_ electrons and the Ny nuclei in the molecule. V(x, ,R,) rep- 
e 


resents the total potential energy of the molecular system under consider- 


ation. The notation used is illustrated in Figure 4.10. 
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Electron 2 Poe a 


Origin 


Figure 4.10 A symbolic sketch of a molecule showing electronic and 


nuclear coordinates 


The potential energy term, V(r,.R), used in the derivation below 


will consist only of the Coulomb interaction terms We? V and ae to 


en 
be defined later. Strictly speaking, other terms like the one due 

to Spin-Orbit interaction, ought to be included if vy) is to represent 
the total potential energy function for the system. However, such terms 
are small (Spin-Orbit coupling, for example, is a relativistic effect) and 
could, if necessary, be treated by Perturbation Theory. It should be noted, 
however, that the Born-Oppenheimer theory is valid regardless of whether 


these terms are neglected or not [78]. 


We write V(r, .R,) as 
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Viena +N ety (4.43) 


where 
Mie is the Coulomb potential term for electron-electron interactions, 


V is the Coulomb potential term for electron-nucleus interactions 


en 
and 
e is the Coulomb potential term for nucleus-nucleus interactions. 
Therefore, 
ee (Ty) re Y — = (4.44a) 
(ij) [rp - x, | 
- 
Ae Zh Se 
= = eat -44b 
WENGER) oe = El (4.44b) 
i 0. 
ee tae ee 
AROS =a) eee (4.44c) 
(a8) |R 8, | 
where 
the notations (ij), (io) and (a8) on the summations imply that each 
pair of particles has to be included just once, 
Ls Ze represent the charges on the nuclei a and 8 respectively, 
7 
and 


e is the electronic charge. 


If the Kinetic Energy Operators in equation (4.42) are denoted by 


TA (rj) and Ty)» i.e. 
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and 
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har (4.45b) 


then the total molecular Hamiltonian becomes 
- =, te a x SUPT pees 
CLE Rs) oT) + T&) + Neato) + Van Ty2Ry) 
+ Van y? (4.46) 


Dropping the coordinate dependence of each term, for convenience, 


equation (4.41) becomes, after substitution of (4.46), 


ay + a + vo + ae + va) Ve =E Ve (4.47) 
Defining 

Ho =T +V (4.48a) 

, e ee 

Ho = T +V (4.48b) 
n nn 

He = Vv (4.48c) 
en 


equation (4.47) can be written as 


(Ho +H +H) ¥, = Ey (4.49) 
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The Born-Oppenheimer Approximation can now be stated mathematically. 
The total wave function Ve (equation (4.41)) of a molecular system 


can be approximated, under certain restrictions, to the product 
GR.) = 0° Gk) - VR) 4.50) 
VE a? Q U) eae a uy) a (4. ) 


where 


v" 38.) is the solution to the equation: 
CURE a, 8) Ue ek ay, (4.51) 


and 


v®) is the solution to the equation: 
@~ nes enon 
(I, + Nps +E)wW =Evy (4.52) 


Furthermore, the total energy of the molecular system can, under the 


same restrictions, be approximated as 
E=E- (4.53) 


The energy operators in equation (4.51) and (4.52) are defined as 
the electronic and nuclear Hamiltonians, respectively. 

It should be noted that equations (4.52) and (4.53) give validity 
to the potential energy diagrams for any molecule. Equation (4.52) states 
that the nuclear motion in a molecule takes place in an effective potential 


energy function Ce a En This is the potential function usually drawn 
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in energy level diagrams. The effect of different electronic states is 
incorporated in this effective potential function by the term E° which 

is the solution to the electronic equation (4.51). Renee: equation 
(4.53) states that the total energy of the ee system is the 
eigenvalue of equation (4.52). 

Should the Born-Oppenheimer Approximation break down, then equations 
(4.50), (4.51), (4.52) and (4.53) no longer hold. In that case, the 
potential energy diagram concept, strictly speaking, has no validity either. 

The conditions under which the Born-Oppenheimer Approximation holds 
will now be derived. It will be assumed first that the approximation 
holds (and that, therefore, the nuclear and electronic motions can be 
separated). An analysis will be made to then derive the conditions which 
have to be satisfied in order that the assumption of separable wavefunctions 
be valid. 


Using equations (4.48) we rewrite (4.51) and (4.52) as, 


e en, @s> =. . es e- o> 
(H +H ) p (r, 5k) = f— (RY) wv (r,5R.) 4,54) 
and 
(+E) y@,) = EVR) (4.55) 
respectively. 


. We expand Ve [in (4.41)] as @ series: 


vp = ve ty + ase ee (4.56) 
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Similarly E [in equation 4.41] is written as: 
cA en (4.57) 


We have then to find the conditions under which the second and following 
terms in (4.56) and (4.57) are negligible compared to the first terms 
in these two equations. 

Substituting (4.56) and (4.57) into the Schroedinger equation 


(4.41) yields 


aqey® +p) + yp 40.) = tae 4 p@ 4 ee 
eae (4.58) 
Therefore, 
yey? + Hy 1) sg ig ST ft we pit), Ce), @) ct ety) 


es 6 (4.59) 


Using H = ie +H +H” we find 


n 


yoy" = Eyoy" + [Tv o]v (4.60) 


where 


[T vod represents the commutator operator defined by 
n 


[Tee eset ever 
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Note that equations (4.54) and (4.55) have been used to derive 
equation (4.60). 
Using (4.59) in (4.60) and ignoring terms of second and higher 


order yields, 


[tr vol™ + Hy) = Ey 4 gDyeyn (4.61) 


(1) a) 


(x, .R,) can be written [78] as the sum 


' R) (4.62) 


where the subscripts ¢€,v,e' and v' are quantum numbers introduced to 
distinguish wavefunctions corresponding to different states. 


Using this expansion in equation (4.60) we obtain 


ae) 
yt geaY 


(= n 
a L [T Vor] Vonye 


(4.63) 


Ignoring the second term in (4.63) (since it is of second order) 


and substituting into (4.61) yields 


qi) n e on 
[Tv] Vey Bes Satay? Eat Vane eat (4.64) 
ae (area n (1).e:.n 
" Lae ee Fo ees nt Ve vey 
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* * 


Multiplying both sides of this equation by Vers or , (where 


ie 
the asterisks denote complex conjugation) and integrating over all 


electronic and nuclear coordinates (tT, and ae respectively) results 


in the following equation [78]: 


* * 

e n e. n iL n 
f was Very lT sbi ly dt odt. a ee Eorytt 

my Loe ch) (1) 
Ey 2 tats + Ee Se eat vy! (4.65) 
where 
GO yes 6.14 ave delta—Dirac, functions. 
Ee vv 


Finally, for the case when ec # e'' and v # v'' equation (4.65) yields 


& 7 
e n e n 
fu ryt 9, Maree. oW ly dtidr 
) Be Cie eee ay fn Coen eee out (4.66) 
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eV wt 


Vv 


This is the result being sought since it provides a condition 


for the validity of the Born-Oppenheimer Approximation. If all the co- 


(1) 
€ 


efficients a rye Were small compared to unity Ce << 1) then the. 


first order correction to the wavefunction, nee in equation (4.62) would 
be small compared to the zeroth order term yoy in equation (4.56). In 


that case 
ve =ypy (4.67) 
which is the Born-Oppenheimer Approximation. 


In general, the commutator in equation (4.66) is not zero. It is 


seen therefore that aaa would be large near a region where electronic 
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degeneracy exists. (i.e. where Ey = Eee o) In other words, it can 
be stated that in general near a region of potential energy curve crossing 
the Born-Oppenheimer Approximation breaks down. 

Now, the experimentally verified fact that efficient E-V transfer 
occurs if the potential energy curves of the final and initial states of 
the system intersect, can be understood. It has been shown that near the 
region of crossing the Born-Oppenheimer Approximation fails resulting in 
a mixing of electronic and vibrational wavefunctions. Under these 
circumstances, it is logical to assume that efficient E-V transfer would 
occur. Furthermore, viewed from this light, the fact (again experimentally 
verified) that a loosely bonded collision eines pean the region of 
curve crossing greatly enhances the efficiency of the process, is under- 
standable also. A collision complex having a lifetime of a few vibrational 


periods would allow enough time for the mixing of wavefunctions, predicted 


by the theory, thus increasing the efficiency of the E-V process. 
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CHAPTER V 


CONCLUSION AND OUTLOOK 


The problem of obtaining a population inversion in the vibrational- 
rotational energy aves of the Hydrogen a Teesre is perhaps best suited 
to optical pumping. 

The small fiecociat ton energy of the molecule and the strong electro- 
negative character of the Hydrogen atom make pumping by the electric- 
discharge or e-beam technique difficult. This is all the more so at high 
pressures. The development of high energy-high power Hy infrared lasers 
must await the overcoming of these difficulties. 

On the other hand, it seems feasible that optically pumped Hy lasers 
could be developed using technology available at present. Such lasers 
would, of necessity, be low energy output devices. The advantages of a 
high pressure tunable IR laser in the 2um region of the spectrum, however, 
make the effort worthwhile. 

Optical excitation of a metastable electronic state (of an atom or 
a molecule) which then transfers its energy to vibrations of H, seems to 
be the most promising technique, since the efficiency of an E-V transfer 
process can be quite high under suitable circumstances, as illustrated in 
Chapter IV. In particular, if the collision complex model holds for a 
given system then the efficiency of the transfer process would, ne a rule, 
be high. This is because, as Tully's formulation shows, the model is not 
very sensitive to some of the parameters (e.g. the Landau-Zener Geet ticiennye 

There is scope for extension of the theory to incorporate the effect 


of pressure explicitly. Since the RRKM expression reduces to the Activated 


Complex formulation in the high pressure limit, the basic mechanism of E-V 
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transfer does not change. However, it would be very useful to have a 
quantitative formulation for the effect of pressure on the vibrational 
partial cross-sections. The effect of deactivating collisions would be 
important at high pressures and this ate to be considered too. 

But, perhaps, the greatest scope is presented when the whole E-V 
ranet ee process is viewed in terms of the breakdown of the Born-Oppenheimer 
Approximations and consequent eck of electronic and nuclear motions. 

It is known that this occurs near a region of electronic degeneracy. It 

would be fruitful to formulate general criteria for the breakdown of this 
approximation. The criteria would probably be based on symmetry considerations 
It would not be necessary to Aspe completely for either the wavefunctions 

or the Hamiltonian for the interacting system of particles. On this basis, 

and making the further restriction that the system form a loosely bonded 
complex, various possible atoms and molecules might be investigated as 


prospective candidates for pumping various molecules including H (E-V 


2° 
transfer between molecules has been demonstrated for the first time, by 
Slanger and Black, recently [80]). 

To illustrate the above points, consider a two level system. 

It can be shown [81] that the energies of such a system are given by 


iy 2 2,12 
patty es MCE es chy ee 4/4) 51°) (54) 


where 
H = H|v = f Hy oe d 
< Sie v. Hy, an 
one v. | : > 
« is the volume over which the integral is evaluated. 


In this case a necessary, but insufficient condition for degeneracy is 
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< ¥, [1H] ¥y> = 0 (5.2) 


If the symmetries of Vy> os and the Hamiltonian are known then it 
would be possible to see whether or THA for a given system of interacting 
particles, settee (5.2) holds. It is not necessary to have complete 
information on these Weare ee only the symmetry properties have to be 
known. 

A development of these Meas ebearine in mind that there exists a 
close relationship between the breakdown of the Born-Oppenheimer criterion 
and efficient E-V transfer should offer a powerful basis for selecting any 


atom or molecule to pump another by this process. 
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APPENDIX I 


SOME TUNABLE LASERS IN THE INFRARED 


Some of the tunable semiconductor lasers in the IR region of the 


spectrum which are presently available are shown in Figure Al below [3]. 


Pb, ony Te 
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Pb, Ge,Te 


€- - oe 


In,.Ga,.,As,Ini F 
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5 10 20 50 100 microns 


Figure Al Wavelength Ranges Covered by Various 


Semiconductor Lasers [3] 
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gas lasers have also been reported by 


These are lasers operating in the 10um 
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APPENDIX II 


VIBRATIONAL BAND FOR Hy 


The table below lists the frequencies and wavelengths of transitions 


between adjacent vibrational levels of Hy for its ground electronic 
4 
(>) ) state. 
g 
The following approximate equation [18] for the vibrational term 


value has been used: 


2 
G W@W) = -w Xv Ziel 
ASO, wav oo ( ) 
where 
v is the vibrational quantum number ,w St Oe, 
and 
wx & 


W xX 
ome) ee 

the values for WO. and WX, being taken as 4935 an and 117 em + 
respectively. Hence, the energy for the transition between levels 


vy = yeand v' =v +1 is 


Gv!) - Gi(v't) = 6 (w#) - Go(v) = w, - ax, LO Woe 


(A2.2) 


Note the anharmonic defect between the lowest two frequencies is 


Av= 0.7 x 1013 Hz 


- —2 
hAv= AE = 2.9 x 10 4 ev = 342x210 ev 
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Table A2: Vibrational Transition Frequencies for Hy 


Transition 
view! 


Frequency 


Wavelength 
x 1013Hz 
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This compares well with the value, 5 x 10m ev obtained from equation 
GID): 


APPENDIX III 


RADIATION TRANSITION PROBABILITIES [18] 


In the presence of electromagnetic radiation the total Hamiltonian 
for an atom or molecule must include the enerEy term due to interaction 
between the radiation and system of particles. 

For instance, if M is the dipole moment of the atom or molecule 


defined by 


where 


ae, laeeand a_ are the unit vectors in the sf and z directions 
respectively,then. the dipole moment interaction 
energy term is M-E (E being the electric field 
intensity vector of the radiation). 

The matrix elements of the electric dipole moment for energy levels 


m and n are defined as 


as . (A3.2a) 
= e a 

R fu, Mead | 

mn * 
= = d (A3.2b) 
é f va MY T 
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uot (A3.2c) 
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where 
ve? ve are the eigenfunctions of the energy levels under 
consideration 

and 


the integration is carried over all space. 


The matrix element 


Rove =cacRe) Flasks aoe 
SOexX yy PA 4 


defines the dipole transition probability for a mn or nm transition. 

It can be shown that if R”” is zero then the transition is forbidden 
as a dipole transition, otherwise not. Magnetic dipole, quadrupole and 
tnduced dtpole transition probabilities are defined in a similar manner. 
In these cases, the corresponding quantity - magnetic dipole moment, 
quadrupole moment or induced dipole moment (due to either the efi or cid 
Berecrs or a combination of the two) is substituted in equations (A3.2) 
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APPENDIX IV 


MARX GENERATOR CIRCUIT 


A ten stage Marx generator, shown in Figure A4.1, was designed 
for the electric discharge experiments described in Chapter II. 


To prevent arcing, the whole circuit was put in a tank of transformer 


oil. 


250K R R R OUTPUT 


re A4.1 Ten Stage Marx Generator 
100 


Figu 
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The output voltave from the “lars as measured across a Fbo dilutes 
Copper Sulphate resistor fs shown ino Uivure Awe, 

Control over the output voltase pulse maynitude was obtained by 
varying the pressure of Ny in the spark gaps. For this best results (an 


output range of between 100 and 179 KV) were obtained when the sides of 


spark gaps containing the triggeriny pins were erounded. 


TOP: 
Charging Pulse 
10 Kv/em 


BOT TOM: 
Output Current 
100A/em 


100 ns/cm 


Figure A4.2 Charying and Output Voltage Pulses trom the Marx 


For the discharge experiments described in Chapter II the Marx 


was fired at a rate of about one pulse per second or slower. 
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APPENDIX V 


CHEMICAL KINETIC NOTATION [57] 


The following list defines the more important of the symbols used in 
Chapter IV. An attempt has been made to correlate the present notation 
with that of other authors, the latter being given in brackets following 


our notation. The authors’ names are abbreviated as follows: 


EiotKk.d, .Laidlerser @lnto0] 
M->R.A. Marcus et al [59] 


teed 6G. VLOLLy [60] 


QQ etc Partition functions [L: Q or q; M:P] 


Ss Number of vibrational degrees of freedom of a molecule 
v Vibration frequencies 

Vv Vibrational quantum numbers 

1 Moments of inertia 

J Rotational quantum numbers 

E Energy per molecule [L: e€ or E; T:¢] 

E. Zero-point energy of a molecule [M:E ;Lie or EO] 

P(E.) Number of quantum states of a given system at the 


quantized energy level E 


k Rate constants in general 

kee First-order rate-constant [L:k'] 

uni 

k Limiting high-pressure value of Saag (M and L: k ] 
Second-order rate-constant [M:k ] 

Sin 5 fe) 


102 


n& hades shane — . 
rot seinen. ated elraie 
aleuliwl eases eegey: 
tial 9, ea 


‘s wy famtr poet eae é apes sage | 
wit ‘s mies ane & x jae " r 


a eae ode er ee oe 
; ib) ah i, a \, - ve 
| "a ea Soa dh SB ait leg aca Pr ae 
adda ; a ye seeds - 5 . | ra 
red oD eae he ) i 


k_ (E*) 


Arr 


103 


Rate-constants (second-order) for energization and de- 


energization, respectively [L:k,,k_ 


Re 1! 


Rate-constant for conversion of energized molecules 


to activated complexes in RRKM theory [Mik 3 L:k, or kp 


+ 
rek(e,J,5)! 
Critical energy for reaction; ‘activation energy’ 
[M:E 5 Lie, or EO] 


Arrhenius activation energy 


A EE Limiting high-pressure values of A and E : and low- 


Py /2- 
* 

hae 
%& 

E 


Arr 


pressure value of Ear 


Pressure at which k = i, 
uni 2 


ry 


Energized molecule and activated complex (L:A” A 
Total non-fixed energy in the active degrees of freedom 
of a given energized molecule AY Lieve, T:W ] 

Total non-fixed energy in the active degrees of freedom 
of a given activated complex a (EY = E. + E. + ET) 
ease T:E* + Ww] 

Vibrational, rotational and translational (along the 


+ + 
reaction coordinate) components of E [for E.3 


= # 4 
Mexs or EES Lie. ; T Ea 
+ 
a = ES + oa (M: Ens nee 
vr Vv r n n 


Energy of adiabatic rotations in the Jth energy level 
* 4: é 
in A and A _ respectively 
| Z 
AEM = Eee Ei Meee Ets Lise ..Ger, 
[ a hi > es oe 
Density of quantum states of A (number of energized 


* k kok 
molecules per unit energy) at energy E [T:N ; M:N (E ) 
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a : + 
) Density of quantum states of A having energy cl 


, ; a 
in the active degrees of freedom and energy Ee along 


the reaction coordinate [M:N (E'-x); LIN, (e7- 7) 


Density of quantum states for the translational motion 


Irn ; : : + 
of A in the reaction coordinate with energy E. 


*)] 


[L:N, (Ce, 


Number of vibrational-rotational quantum states of re 
at the quantized energy level E [M:P (") or a (ET); 
T:N'] 

The arbitrary length of region at the top of potential 
energy barrier which is taken to define the activated 
complex [M:b] 

Characteristic mass for motion in the reaction 
coordinate [M:msL:M" ] 


; “+ 
Rate constant at which complexes of energy E. cross the 


barrier [Mik,3 L3x] 


; ; oF 
Partition functions for addiabatic rotations in A and A 


: pte ds : 
respectively [M:P,,P, 5 L:Q,> Qe ] 


Partition functions for active degrees of freedom A 


and An respectively (at has one fewer degree of freedom) 


ee ed # 
[M:P,, Pos LiQ,, Q, 


Complete vibrational-rotational partition functions for 
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APPENDIX VI 


A NOTE ON PARTITION FUNCTIONS [57, 58] 


Definition 
For a molecular system capable of existing in a series of 


quantized energy levels with energies E Ene see Che COtal partition 


oe Ae 


function for the molecule is defined as 


CO 
i} 


g exp (-E)/KT) +8, exp (-E, /KT) 1» Sp, exp (-E,/KT) Sera 


2 8, exp(-E,/kT) (A6.1) 
i=o 
where 
5 is the degeneracy (or statistical weight) of the energy 
level E> which may be defined as the number of different 
independent wave functions of the system with energy E, 
(i.e. the number of physically distinct quantum states of 
that energy). 


If the vibrational, rotational and electronic motions in a molecule 


can be separated then 


=E + +E Fee 
ED =A ey r t 


where 


Ee is the total energy of the molecule 


Een be, jes. and E. are the electronic, vibrational, rotational 
e r 


and translational energies respectively. 
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Under these conditions it can be shown [58] that the total partition 
function, Q, may be written as a product of the vibrational, rotational, 


electronic and translational partition functions. 


Cea: Oa, (A6.2) 


Note that since nuclear and electronic motions are only separable under 
the conditions of the Born-Oppenheimer approximation,equation (A6.2) 


holds only if the approximation does. 
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APPENDIX VII 


COMPUTER PROGRAM FOR CALCULATING VIBRATIONAL CROSS-SECTIONS 


The Rhomberg method [70] of integration was used to evaluate the 
total quenching and partial vibrational cross-sections for the de- 
excitation of Na" by collisions with ei 

Figure A7.1 shows the flow chart for the complete program while 


Figure A7.2 shows the flow chart for the Rhomberg sub-routine. 
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Figure A7.2 Block Diagram of Program Used 
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Define flx) 
Read a,b,€,0 


X,=a 
h=b-a 


T, =(h/2) Ifla) + f(b)] 


Figure A7.2 The Rhomberg Subroutine [70] 
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APPENDIX VIII 


A NOTE ON ATOMIC UNITS 


Since atomic quantities like electronic ae mass and energies 
are very small when expressed in Sudtnary waite (Bore oe orm, kes) sit5is 
Hecesatane for convenience, to define a system of “Atomic Units". There 
are many choices available such that the values of fundamental constants, 
as expressed by these units, are conveniently - handled small numbers. 
Any system could be used, provided it is based upon the least number of 
independent quantities [82]. 

For the computations involved in Chapter IV, the following basic 
units were defined. 

Mass: (1)atomic mass unit (a.m.u.) 

Energy: (1)electron volt (ev) 

i ° ° 

Length: (1) Angstrom Unit (A) 

The units of the other quantities required for computations were 
defined by using the above three units. 

Hence, the units of van der Waals' constant, angular momentum and 


cross-sections used were, 


van der Waal's constant: (nie 
Angular momentum: (eyyt!2 (asmeusyt!* (a) 
Cross-section: (ay? 

The following conversion factors apply: 

1 a.m.u. = erlOmes bam 

1 ev = 1260 x ane erg 

ib (ev) (A) = 1.60 x TORGe Gre) (ae 
LA = 107° em 

(ev) )/2 (a.m.u.)1/2 dy = 1 Gu We sec 
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